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This  paper  describes  a  model  for  survival  functions  called  the  cumulative 
damage  threshold  crossing  model*  Under  this  model,  an  item,  consisting  of 
a  large  number  of  components  which  suffer  damage  at  regular  moments  of  time, 
fails  as  soon  as  the  maximum  cumulative  damage  to  some  component  crosses  a 
certain  threshold.  Under  various  conditions  on  the  damage  distribution,  two 
plausible  new  survival  distributions  , 

Sj(t)  »  exp{-exp(a+Bt,s-Yt~,5)}»  t  >  0,  -“  <  a  <  0  >  0  and  y  >  0, 


and 


S2(t)  ■  exp{-Bt(l/t  -  Y)”a)»  0  <  t  <  1/y»  o  >  -1,  &  >  0  and  y  >  0 

0*9  (  \ 

arise.  -She'  survival  function  CO  is  shown  to  fit  several  data  sets  in  the 
area  of  cancer  studies,  failure  times  of  insulating  fluids  and  failure  times 
of  air  conditioning  units. 


1.  Introduction. 

The  achievements  of  Reliability  Theory  depend  largely  on  successful 
modelling  of  survival  functions  of  lifetimes.  The  famous  'shock  model*  of 
Esary,  Marshall  and  Proschan  (1973)  can  be  described  as  follows.  An  item 
is  buffeted  by  shocks  which  arise  according  to  a  Poisson  process.  Each  shock 
causes  a  random  damage  and  the  damages  accumulate.  The  item  fails  as  soon  as 
the  cumulated  damage  exceeds  a  certain  threshold.  Further  generalizations 
of  this  model  have  appeared  in  the  literature,  e.g.  A-Hameed  and  Proschan 
(1973,  1975).  Only  qualitative  properties  of  survival  functions  under  this 
model  were  studied  in  these  papers.  The  survival  functions  were  not  fitted 
to  actual  data. 

The  classical  results  of  Gnedenko  (1943)  on  the  three  types  of  asymptotic 
distributions  for  an  extreme  observation  lend  themselves  naturally  to  the 
modelling  of  survival  functions.  That  is  why  in  many  applications,  survival 
functions  are  assumed  to  be  Weibull 

(!•!)  S4(t)  ■  exp{-(uc)*},  t  >  0,  a  >  0,  0  >  0, 

r 

or  to  be  the  shifted  Weibull 

(1.2)  S3(t)  -  exp(-[a(t-y))8},  t  >  y,  a  >  0,  0  >  0,  -•  <  y  <  «. 

There  are  many  variations  to  the  use  of  Gnedenko's  (1943)  results  in 
modelling  survival  distributions.  One  can  introduce  dependence  among  the 
observations,  one  can  have  spare  parts  which  instantaneously  replace  parts 
that  fail  and  thus  increase  system  life,  etc.  A  partial  list  of  references 
that  contain  such  extensions  is,  Ashar  (1960),  Bimbaum,  Esary  and  Saunders 
(1961),  Davis  (1952),  Flehinger  (1962),  Harris  (1970),  Laurent  (1958), 

Schafer  and  Finkelstein  (1962),  Sen,  Bhattacharyya  and  Suh  (1973)  and  Weiss 
(1961) . 


In  this  paper  we  suggest  a  different  aodel  which  can  be  described  in 
broad  tens  as  follows.  An  iten  consists  of  a  large  number  of  components. 

At  regular  moments  of  time,  t  ■  1,  2,  . ..,  random  damages  occur  to  the  com¬ 
ponents  and  the  damages  accumulate.  The  item  fails  as  soon  as  the  cumulative 
damage  to  some  component  crosses  a  certain  threshold. 

It  is  easy  to  produce  examples  where  such  a  model  would  apply.  In  cancer 
studies  one  may  define  the  onset  of  cancer  as  the  time  at  which  one  of  the 
cells  becomes  mutant.  Thus  one  may  postulate  that  an  organism  consists  of  a 
large  number  of  cells,  that  the  environment  produces  random  malevolent  effects 
on  the  cells  and  that  a  cell  becomes  mutant  when  the  cumulative  malevolent 
effect  exceeds  a  certain  threshold.  Similar  models  can  be  postulated  in  other 
situations  like  machines  with  a  large  number  of  parts,  electronic  equipment 
with  large  numbers  of  transistors,  diodes,  etc. 

For  the  simplest  form  of  our  model,  we  present  some  asymptotic  results 
in  Section  2.  This  leads  to  two  new  survival  functions 

(1.3)  Sj (t)  ■  exp(-exp(  o+Bt^-Yt"15  )},  t  >  0,  -•<»<•,  p  >  0,  y  >  0 

and 

(1.4)  S2(t)  ■  exp{-Bt(l/t-  y)"°),  0  <  t  <  1/y,  o  >  -1,  B  >  0,  y  >  0. 

the  proofs  of  these  results  employ  results  on  rates  of  convergence  for  moderate 
and  large  deviations. 

In  Section  3  we  fit  our  survival  functions  Sj(0  and  S2(»)  as  well  as 
the  classical  ones  Sj(*)  and  S^(>)  to  a  variety  of  data.  The  data  come  from 
cancer  studies  in  England  and  Wales,  Doll  (1971),  from  studies  on  insulating 
fluids,  W.  Nelson  (1972,  197S),  and  from  air  conditioning  systems  on  Boeing 
aircraft,  Proschan  (1963).  It  is  gratifying  to  see  that  the  new  survival 
function  S^(*)  fits  nearly  all  these  data  better  than  the  others. 


The  survival  functions  Sj(»),  S2(»)  and  S^C*)  each  contain  three  parameters 
a,  0  and  y  while  S4(*)  contains  only  two  parameters  o  and  0.  The  survival 
functions  S2(«)»  S3(*)  and  S4(»)  may  have  either  an  increasing  failure  rate 
(IFR)  or  a  decreasing  failure  rate  (DFR) .  For  example,  S2(«)  is  IFR  if  o  a  0 
and  Sj(-)  and  S4(»)  are  each  IFR  if  0  a  1.  On  the  other  hand,  Sj(*)  is  a 
more  flexible  survival  function.  Sj(*)  may,  depending  upon  the  values  of 
its  parameters,  be  an  IFR  survival  distribution  or  it  may  have  an  increasing 
failure  rate  for  small  and  large  values  of  t  while  having  a  decreasing  failure 
rate  for  moderate  values  of  t. 

2.  The  survival  functions  SjCO  and  S2(*J. 

Consider  an  item  that  consists  of  k  components.  The  i1*1  component  suffers 

t 

damage  Z,^  at  time  t  ■  1,  2,  ...  .  Let  S.^.  »  V  Z, .  be  the  cumulative  damage 

It  XZ  j.j  ij 

to  component  i  at  time  t,  i  *  1,  ...,  k,  and  let  S*  «*  max  S.  be  the  maximum 

Kt  lsisk 

cumulative  damage  suffered  by  any  component.  The  cumulative  damage  threshold 
crossing  model  that  we  propose  in  this  paper  postulates  that  the  life  of  the 
item  is  given  by 

(2.1)  L  *  inf{t:  S*t*d), 

where  d  is  the  threshold. 

We  make  the  following  simplifying  assumption.  The  random  damages 
Zjj,  Zj2»  • . . ,  Z2j ,  Z22,  ...,  ...,  Z^2»  ...  are  independent  and  identically 

distributed  with  common  distribution  function  F(«).  One  could  relax  this  by 
assuming  that  the  vectors  (^.....Z^),  (Z^2,...,Zj(2),  ...  are  independent 
and  identically  distributed  or  follow  a  special  stochastic  process.  One  could 
also  define  the  life  of  the  item  by 


L  •  inf{t:  ♦(slt , .  -  -  ,Sfct)  id). 


-  4  - 


where  ♦(•,...,•)  is  a  function  of  the  individual  damages  representing  the 
damage  to  the  item.  We  propose  to  examine  these  generalizations  in  a  future 
paper. 

Reverting  our  attention  to  the  simple  model,  we  obtain  the  asymptotic 
distribution  of  as  t  +  »  and  k  ■  k(t)  +  »  at  a  certain  rate  in  Theorems  1 
and  2.  In  Theorem  1  we  assume,  among  other  conditions,  that  F(»)  has  a  finite 
variance  and  obtain  a  double  exponential  asymptotic  distribution.  In  Theorem  2 
we  assume,  among  other  conditions,  that  F(*)  has  a  support  that  is  unbounded 
above  and  a  tail  that  decreases  to  0  like  a  power  of  x  *  and  obtain  a  Weibull 
asymptotic  distribution.  By  inverting  these  distributions  to  obtain  the 
approximate  distribution  of  the  lifelength  L  we  are  led  to  the  survival 
functions  Sj(»)  and  S2(*)  of  (1.1)  and  (1.2). 

Theorem  1.  Let  the  mean  y  and  variance  o2  of  F(*)  be  finite.  Further,  for 
some  q  in  (0,1),  let  F(0  satisfy 

(2.2)  1  -  F(z)  -  5 (log  z)(1+q  )/2) 

and 

<0 

(2.3)  /  u2dF(u)  *  o(l/log  z) 

z 


as  z  +  ®.  Let  k  =  k(t)  be  an  integer  valued  function  satisfying 


(2.4) 


kst 


q2/z 


for  large  t.  Set 

(2.5)  «kt  ■  o(t/2  log  k)1* 


and 


» 


» 
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(2.6)  ^  •  yt  ♦  o(2t  log  k)**  -  (l/2)«»kt (log  log  k  ♦  log  4ir). 

Then,  for  -•  «  y  <  •,  u  t  +  », 

(2.7)  P(S*t  s  *  exp(-exp(-y)}. 

Proof:  Since  the  ' s  have  finite  variance,  the  random  variables  S^t  are 
approximately  normally  distributed.  The  maximum  of  normally  distributed 
random  variables  has  an  asymptotic  double  exponential  distribution.  Condi¬ 
tions  (2. 2) -(2. 6)  allow  us  to  approximate  moderate  deviation  probabilities 
of  S^t  adequately  and  then  to  use  standard  extreme  value  theory  to  establish 

(2.7) ,  as  shown  below. 

Theorem  1  of  Amosova  (1978)  on  moderate  deviations  states  that  under 
conditions  (2.2)  and  (2.3), 

(2.8)  P{SU  >  ut  +  xot*}  -  (l-*(x))(l*o(l)) 

uniformly  for  x  in  (0,  q(log  t)**),  as  t  +  ■,  where  *(•)  is  the  standard 
normal  distribution  function. 

Notice  that  under  condition  (2.4),  for  any  y, 

°ytr  ♦  ^  •  lit  . 

where 

(2.9)  xkt  «  (2  log  k)*5*  (2  log  k)"*{y-  (l/2)(log  log  k  ♦  log  4*)>  <  q(log  t)** 
for  large  t.  Thus 

kP{Slt  >  «kty*Bkt>  “  Ml-*(xkt))(l^°(l))  *  e’y 

as  t  >  •  from  (2.8)  and  a  standard  result  on  the  maximum  from  the  normal 
distribution,  for  instance  see  Galaabos  (1978,  Sec.  2.2.3).  Hence 


-  6  • 

1 i 

log  Pfsjt  S  okt1,*skt)  -  k  log  Mslt  S  v>'*0kt) 

'  -kp(sit  ”  V  *■*>  * 

9 

which  implies  (2.7).  □ 

In  the  next  theorem  we  assume  that  the  support  of  F(«)  is  unbounded 

»  ~ 

above  and  has  a  tail  that  goes  to  zero  like  a  power  of  x'1.  A  distribution 

F(*)  can  satisfy  the  conditions  of  both  Theorem  1  and  Theorem  2.  For  such 

F(*)»  one  can  obtain  two  different  asymptotic  distributions. because  k  «  k(t) 

may  tend  to  •  at  different  rates. 

» 

Theorem  2. 

Let  F(«)  have  a  finite  mean  u  and  satisfy 

»  • 

(2.10) 

1  -  F(x)  =  x”°h(x) 

where  h(*) 

is  a  slowly  varying  function  and  a  >  1.  Let 

»  < 

(2.11) 

Yt  *  inf{x:  x’ah(x)  s  1/t) 

and 

>  j 

(2.12) 

«t  *  vt. 

Let  k  «  k(t)  •  so  that 

t  ( 

(2.13) 

ta/k  o  0(1) 

at  t  +  •. 

Then,  for  x  >  0, 

»  .  < 

(2.14) 

MSgt  *  Ykt**«t>  ♦  exp(-x"a) 

as  t 

1  ( 

1  i 

-  7  - 

Proof:  Condition  (2.10)  implies  that  the  support  of  F(*)  is  unbounded  above. 
Hie  fact  that  a  >  1  implies  that  F(«)  has  a  finite  mean,  there  are  two  more 
uses  of  Condition  (2.10)  in  the  proof.  The  asymptotic  distribution  of  a 
maximum  from  F(»)  will  be  exp(-x"a)  (see  for  instance  Galambos  (1978,  Theorem 
2.4.3));  in  fact,  for  x  >  0, 

(2.15)  log  Ft(Yf.x)  -  t  log  F(ytx)  *  -t(l-F(Ytx))  -  -x“a 
as  t  ■»  ». 

The  second  use  of  (2.10)  comes  in  approximating  large  deviation  proba¬ 
bilities  of  Slt .  Let  q  >/V  Theorem  1  of  S.  V.  Nagaev  (1982)  states  that 
under  Condition  (2.10) 

(2.16)  P(Slt>x)  «  t(l-F(«))(l*o(l)) 

as  t  +•,  uniformly  for  x  >  qt.  From  the  definition  of  yt  in  (2.11)  and  from 
(2.15)  it  can  be  seen  that 

(2.17)  Yt  -*■  ®  as  t  ®  and  Yt  =  t^as(t) 
where  s(t)  is  slowly  varying.  Thus 

(2.18)  Ykt  -v  k1/aYt  . 

Condition  (2.13)  relating  k  to  t  can  be  rewritten  as  k^a  2  pt  for  some 
p  >  0.  Hence,  for  any  x  >  0  and  q  >  0, 

Yktx  ♦  6t  *  kl/0Yt*  ♦  t|i  >  PtYtx  ♦  Ut  >  qt 

for  large  t  since  Yt  •.  From  (2.16)  and  (2.15),  it  follows  that,  for 
x  >  0,  log  P(S*t  s  Yktx>«t>  ■  k  log  P(Slt  s  Tktx*6t>  -v  -kP{Slt  >  Yktx*«t} 

*  -kt(l-F(Yktx))  *  -x*a  as  t  ♦*.  This  establishes  (2.14).  0 


» 


9 


» 


I 
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We  will  now  show  how  the  survival  functions  S^C*)  and  S2(*)  of  (1.1) 
and  (1.2)  arise  naturally  from  Theorems  1  and  2.  From  the  definition  of  L 
in  (2.1), 

(2.19)  P(LStt)  =  P{S*t<d). 

Under  the  conditions  of  Theorem  1,  this  is  approximated  by 
exp(-exp((d/o)  ((2  log  lO/t)*5-  (y/o)  (2t  log  k^+gOO)) 

where 

g(k)  «  2  log  k  -  (1/2) (log  log  k  ♦  log  4*). 

2 

Since  log  k  s  (1/2) q  log  t,  we  can,  after  ignoring  terms  involving  log  t, 
approximate  the  survival  function  by 

exp{-exp(a  ♦  gt*5  -  yt**)  } . 

The  parameters  a,  $,  y  satisfy  -«  <  a  <  »,  g  >  0  and  y  »  0,  since  in  practice 
one  would  assume  that  ECZ^)  >  0  and  d  >  0.  This  is  the  survival  function 

S.  ( •) •  Loosely  speaking,  the  assumptions  required  for  this  are  that  the  damage 

1  2 

2  a  /2 

distribution  P  possess  moments  of  order  2  ♦  q  and  that  k  tend  to  «®  like  tH 
or  slower. 

Under  the  conditions  of  Theorem  2,  by  specializing  condition  (2.13)  to 

read 

k  ■  pt° 

for  some  finite  p,  the  survival  function  in  (2.19)  can  be  approximated  by 

d-6  "°  ♦%- 

exp|- (-^7)  }  *  ®xp|-P°^ta  (d-yt)’°/s(t)| 
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where  s(t)  is  a  slowly  varying  function  of  t.  As  before  we  ignore  the  slowly 
varying  function  and  approximate  the  survival  function  by 

exp{-$t(l/t-Y)"a> 

for  0  <  t  <  1/y  and  where  0  >  0,  y  >  0  and  a  >  1.  This  leads  to  the  survival 
function  Loosely  speaking  the  assumptions  required  for  this  are  that 

the  tail  of  the  damage  distribution  be  slowly  varying  of  order  «  >  1  and 
that  k  tend  to  00  like  t®  or  faster. 

We  will  now  give  a  counterexample  to  Theorem  1  to  show  that  when  condi¬ 
tion  (2.4)  relating  k  to  t  is  violated  the  asymptotic  distributi  >f  S*t 
need  not  be  the  double  exponential  or  even  one  of  the  other  two  treme  value 
distributions.  Let  Zll*  Z12*  -  be  independent  and  identical:  i' ^tributed 

with  P(Zjj»0)  ■  PCZjj-l)  «  1/2.  Let  k  >  ^log  2],  where  [•]  ±s  the  largest 
integer  function.  Then  conditions  (2.2)  and  (2.3)  are  fulfilled  but  not  con¬ 
dition  (2.4).  Note  that  S£t  is  integer  valued  and  satisfies  0  s  S*t  s  t. 
Further, 

P{S*tit-l)  -  (1-1/2*)*  1/2  , 

and 

P{S*tst-2)  -  (1-  (t+n/zV  -V  0  , 

as  t  Thus  the  limiting  distribution  of  S*t  -  t  exists  and  takes  on 

values  -1  and  0  with  probabilities  1/2  and  1/2. 

In  the  above,  by  using  Theorem  2,  we  have  given  a  model  for  the  survival 
function  S2(»)  only  when  a  >  1.  However,  it  is  clear  from  the  form  of  S2(*) 
that  it  is  a  survival  function  for  a  >  -1.  When  we  fit  this  distribution 
in  Section  3  we  will  allow  the  range  of  a  to  be  (-!,<»). 


10  - 


By  computing  derivatives  it  is  easy  to  see  that  S2(*)  is  IFR  for  a  >  0, 
exponential  for  a  >  0  and  DFR  for  -1  <  a  <  0.  Similarly,  it  can  be  shown 
that  SjC*)  is  IFR  when  By  *  (9*6/T)/16  -  1.212.  If  By  <  (9«-6«/I)/16,  then 
the  failure  rate  of  Sj(»)  initially  increases,  then  decreases  and  finally 
becomes  an  increasing  function. 

3.  Fitting  Sj(*)  and  S^CO  to  data. 

Many  authors  have  proposed  probability  models  to  explain  observed  cancer 
incidence  rates.  Creasy  (1981)  ascribes  the  most  popular  model  to  Armitage 
and  Doll  (1961)  in  which  a  multistage  theory  of  cancer  development  was  pro¬ 
posed.  This  leads  to  the  Neibull  survival  function  S^(>).  Later,  Doll  (1971) 
extended  this  model  to  allow  for  an  initial  dormant  period,  which  lead  to  the 
shifted  Neibull  survival  function  S3(*).  The  cumulative  damage  threshold 
crossing  model  proposed  in  this  paper  suggests  two  new  survival  functions, 
S1(«)  and  S2(*)- 

Doll  (1971)  gives  data  on  incidence  rates  (x  100,000)  in  age  groups  of 
0-4,  S-9,  ...,  70-74  years,  for  several  types  of  cancer  in  England  and  Wales 
from  1961  to  1963.  Denoting  the  failure  rate  of  5^(0  by  r ^ ( •) »  we  estimated 
the  parameters  of  S^(«)  by  minimizing 

16  , 

(3.1)  l  (r(t.)  -  r. (t.)j 

i=l  1  J  1 

where  t^  is  the  midpoint  of  the  i^  age  group  and  r(t^)  is  the  corresponding 
observed  incidence  rate.  The  Nelder-Mead  algorithm  was  used  for  this  mini¬ 
misation,  as  in  L.  S.  Nelson  (1973),  which  included  the  constraints  on  the 
parameters.  The  residual  sum  of  squares,  namely  the  minimum  value  attained 
by  (3.1),  will  be  denoted  by  R^,  j  =  1,  ...,  4.  This  is  to  be  compared  with 


-li¬ 


the  corrected  sua  of  squares,  Rq,  given  by 


(3.2) 


l  (r(tj  -  r) 
i-1  1 


where  r  is  the  average  observed  incidence  rate. 

The  data  were  classified  into  seven  types  of  cancer  as  indicated  in 
Table  1.  The  category  'other'  represents  cancer  of  the  bone,  testes,  rectus 
or  prostate.  Table  1  gives  a  summary  of  these  results.  For  all  cancer  types 
except  leukemia.  Model  1  provided  the  best  fit. 

N.  Nelson  (1972,  1975)  examined  data  on  time  until  failure  to  insulate 
when  charges  of  several  thousand  volts  were  applied  to  parallel  plates 
separated  by  an  insulating  fluid.  The  parameters  of  the  survival  function 
Sj(*)  were  estimated  by  minimizing  the  von-Mises  statistic 


?  <pn<V-W>2 

itl  ijCt.Hl-s.cV) 


where  tj,  ...»  tft  are  the  observed  times  of  failures  arranged  in  increasing 
order,  and  Pn(t.)  *  .  The  minimum  value  of  (3.3)  will  be  denoted  by 

Vj»  j  *  1»  ...»  4.  The  results  are  reported  in  Table  2.  Once  again,  Sj(») 
seems  to  give  the  best  fit. 

The  third  set  of  data  analyzed  comes  from  Proschan  (1963)  on  times  to 
failure  of  Boeing  720  aircraft  air  conditioning  units.  The  following  modified 
von-Mises  statistic  was  minimized  in  order  to  estimate  the  parameters  of  S ^ ( •) 
and  S3(«): 

n  cyy-yy)2 

i-l  G1|Cti)811CtJ) 


(3.4) 


where  again  tj,  ....  t  are  the  observed  tines  of  failure  arranged  in  in¬ 
creasing  order,  and  fi  (t.)  *  1  -  G  (t.)  ■  -----  .  The  minimum  value  of  (3.4) 

n  i  n  l  n 

will  be  denoted  by  W^,  j  «  1,  3.  The  results  are  reported  in  Table  3  and 
Sj(*)  gave  the  better  fit. 

The  first  data  set  on  cancer  appears  to  have  an  IFR  distribution,  while 
the  second  and  third  data  sets  indicate  DFR  distributions.  In  fact,  Proschan 
(1963)  rejects  the  hypothesis  of  constant  failure  rate  in  favor  of  a  DFR 
distribution  for  the  pooled  failure  times  of  all  of  the  air  conditioners. 

As  fitted  to  the  air  conditioner  data,  Sj(t),  has  a  decreasing  failure  rate 
on  the  observed  range  from  t  -  11.5  to  t  ■  231.7.  It  is  heartening  to  see 
that  Sj(*)  fares  well  in  all  the  cases  considered. 
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Table  1.  Results  of  fitting  S.(«)»* ••»$.(•)  to  cancer  incidence 
(1971) . 


Sj(-) 

parameter 
estimates 


data  from  Doll 


1 

s4(0 

Ri  parameter  1 

R*  parameter 

estimates  | 

estimates 

Leukemia  (1862  males ,  R*  *  1.61E-7) 


5.21E-9 

-l.S2E*l 

1.14 

1.00E-4 

6.38E-9 

4.66E-9 

4.14  2.95E-3 
1.00E-12  1.52E+1 
1.01E-4  J  -1.64E+2 

6.39E-9 

4.69E-3 

5.16 

Stomach 

(9274  males,  R0  =  8.34E-6) 

3.80E-8 

6.80 

6.30E-6 

8.93E+1 

8.27E-8 

|  7.86E-8 

4.0S  I  6.79E-3 

1.00E-11  j  4.93 

7.57E-5  |  2.50 

8.22E-8 

6.71E-3 

5.09 

Skin  (9012  males,  RQ  »  6.89E-6) 

8.54E-9 

-6.89 

7.81E-1 

2.98E+1 

1.00E-8 

9.48E-9 

4.09  6.32E-3 

7.61E-12  5.83 

9.14E-5  -1.07E+1 

1.01E-8 

6.64E-3 

5.12 

Lung  (26676  males,  Rq  a  4.88E-5) 

4.60E-6 

5.15 

7.23E-6 

6.55E+1 

5.85E-6 

J  5.76E-6 

2.75  |  7.13E-3 

8.32E-9  J  3.63 

8.35E-8  |  2.50 

5.85E-6 

7.04E-3 

3.76 

Other  (10565  males,  Rq  =  2.25E-5) 

2.41E-8 

1.18E+1 

3.57E-2 

1 . 35E+2 

2.97E-8 

[  2.86E-8 

6.75  8.77E-3 

9.30E-17  7.52 

1.04E-4  2.50 

2.93E-8 

8.60E-3 

7.76 

All  Male  (57389  males,  RQ  =  2.89E-4) 


2.34E-6 


4.41E-6 


4.23E-6 


4.39E-6 


7.98 

4.82E-7 

8.36E+1 


3.73 

2.56E-10 

8.94E-5 


9.62E-3 

4.61 

2.50 


9.41E-3 

4.76 


Breast  (21017  females,  Rq  ■ 

7.76E-6) 

1.89E-7 

3.96E-7  I  2.24E-7 

3.96E-7 

2.53 

1.83  1 

4.35E-3 

4.66E-3 

4.49E-10 

2.47E-7  | 

2.16 

2.83 

4.77E+1 

5.79E-9  1 

1.75E+1 

-- 

*RQ  is  the  corrected  sum  of  squares,  (3.2),  and  R. ,  j  •  1,...,4  are  the 
minimized  residual  sum  of  squares,  (3.1).  * 
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Table  2.  Results  of  fitting  Sj (•)»•• • ,S^(*)  to  data  on  failure  times  of  insulators 
from  W.  Nelson  (1972,  1975). 


s. 

*5 

L(0 

parameter 

estimates 

s. 

4 

(• ) 

parameter 

estimates 

S2 

*3 

[(*) 

parameter 

estimates 

S4 

n 

parameter 

estimates 

0.107 

-0.523 

0.140 

0.923 

0.238 

32  KV  ( 

-0.513 

0.235 

1.469E-3 

n  ■  15) 

0.191 

0.051 

0.424 

0.218 

0.205 

0.047 

0.464 

0.341 

-0.317 

0.231 

1.S12 

0.296 

34  KV  (n 

-0.283 

0.190 

6.357E-3 

«  19) 

d.295 

0.095 

0.675 

0.082 

0.303 

0.094 

0.692 

0.200 

1.200 

0.183 

2.640 

0.315 

36  KV  (n 

-0.058 

0.321 

3.738E-2 

*  15) 

0.249 

0.335 

0.836 

0.256 

0.328 

0.295 

0.940 

*Vj ,  j  »  1, . . . ,4  are  the  minimum  values  of  the  von-Mises  statistic,  (3.3). 


» 


I 


Table  3.  Results  of  fitting  Sj(«)  and  Sj(»)  to  213  failures  of  Boeing  air 
conditioners  from  Proschan  (1963). 


Sj(-) 

Si 

><*> 

W? 

Parameter 

K 

Parameter 

1 

estimates 

estimates 

0.3055 

0.2945 

0.5532 

0.1170 

0.0728 

0.8978 

8.7581 

1.0000 

| 


*Wj  and  are  the  minimum  values  of  the  modified  von-Mises  statistic,  (3.4).  • 


I 
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and 
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arise.  The  survival  function  Sj(0  is  shown  to  fit  several  data  sets  in  the”  area 
of  cancer  studies,  failure  times  of  insulating  fluids  and  failure  times  of  air 
conditioning  units. 


